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An approach is described to calculate Generalized Parton Distributions (GPDs) in 
Constituent Quark Models (CQM). The GPDs are obtained from wave functions to be 
evaluated in a given CQM. The general relations linking the twist-two GPDs to the form 
factors and to the leading twist quark densities are recovered. Results for the leading 
-^j ■ twist, unpolarized GPD in the Isgur and Karl model are presented. 

q_{ 1. INTRODUCTION 

<d : 

Generalized Parton Distributions (GPDs) are becoming one of the main topics of inter- 
est in hadronic physics M. GPDs are a natural bridge between exclusive processes, such 
as elastic scattering, described in terms of form factors, and inclusive ones, described in 
terms of structure functions. As it happens for the usual Parton Distributions (PDs), the 
measurement of GPDs provides us with a unique way to access several features of the 
structure of the nucleon, such as the quark orbital angular momentum contribution to 
the proton spin 0|| . Therefore, relevant experimental efforts to measure GPDs will take 
place in the next few years and it becomes urgent to produce predictions for these quanti- 
ties. Several calculations have been already performed by using different approaches 
and an impressive effort has been devoted to study their QCD evolution properties |5|]. 

A step towards calculations of GPDs in Constituent Quark Models (CQM) can be 
found in ||, and a consistent approach has been proposed in 0. The CQM has a long 
story of successful predictions in low energy studies of the structure of the nucleon. In 
the high energy sector, in order to compare model predictions with data, one has to 
evolve, according to QCD, the leading twist component of the physical structure functions 
obtained at the low momentum scale associated with the model. Such a procedure, already 
addressed in ||, has proven successful in describing the gross features of standard PDs by 
using different CQM (see, e.g., [[|). Similar expectations motivated the study of GPDs 
described in J7j, where a simple formalism is described to calculate GPDs from any model. 
In this talk, the approach of [7| is reviewed and applied to the Isgur and Karl (IK) 
model. 



2 



2. GENERAL FORMALISM 

Let us think to diffractive DIS off a nucleon target, with initial and final momenta P 
and P', respectively. GPDs describe the amplitude for finding a quark with momentum 
fraction x + £/2 (in the IMF) in the nucleon and replacing it back into the nucleon with 
a momentum transfer A M . The GPD H(x,^,A 2 ) is introduced by defining the twist-two 
part of the light-cone correlation function |J 

= H(x,Z,A 2 )U(P')rU(P) + E(x,^A 2 )U(P')—-^U(P) + ... (1) 



where ellipses denote higher-twist contributions, ip is a quark field and M is the nucleon 
mass. The £ variable, the so called "skewedness" , is defined by the relation £ = —n ■ A, 
where n = (1, 0, 0, — 1)/(2A) and A depends on the reference frame. The £ variable is 



bounded by and y/-A 2 /^JM 2 - A 2 /A. Besides, one has t = A 2 = A 2 , - A 2 . When the 
longitudinal momentum fraction of the quark is less than — £/2, GPDs describe antiquarks; 
when it is larger than £/2, they describe quarks; when it is between — £/2 and £/2, they 
describe qq pairs. There are two natural limits for H(x,£, A 2 ): i) when P' = P, i.e., 
A 2 = £ = 0, the so called "forward" limit, one recovers the usual PDs 

H(x,0,0) = q(x) (2) 

ii) the integration over x yields the Dirac Form Factor (FF) 

[ dxH(x,£, A 2 ) = Fi(A 2 ) . (3) 

Any model estimate of the GPDs has to respect the above two crucial constraints. 

In Ref. M, the Impulse Approximation (IA) expression for the GPD H(x,£, A 2 ), suit- 
able to perform CQM calculations, has been obtained. In particular it has been found 
that, substituting the quark fields in the left-hand-side of Eq.(l), taking into account the 
quarks degrees of freedom only, using IA, considering a process with A 2 <C M 2 in the 
Nucleon rest frame, using a symmetric wave function (as the one given in a NR quark 
model once color has been taken into account), one obtains 

H(x, £,A 2 ) = fdk 5[x+^~ \- ) n(k, k + A) , (4) 



2 M f 

where n(k, k + A) is the one-body non-diagonal momentum distribution: 
n(k, k + A) = 3 / ^*(ki, k 2 , k + A)ip(ki, k 2 , k)dk x dk 2 = 

= Je i « UK) *-^ ) p{r,r')drdr' , (5) 
defined through the one-body non diagonal charge density 

p(r,r) = / tp*(ri,r 2 ,r)^(ri,r2,r)df 1 df 2 ■ (6) 
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Eq. (4) allows the calculation of the GPD H(x,£,,A 2 ) in any CQM, and it naturally 
verifies the two crucial constraints, Eqs. (2) and (3) @. With respect to Eq. (4), a 
few caveats are necessary. Due to the use of CQM wave functions, only quarks GPDs 
can be evaluated, i.e., only the region x > £/2 can be explored. The approach has 
to be improved in order to study the sea region (— £/2 < x < £/2). Besides, in the 
argument of the 5 function in Eq. (4), due to the used approximations, the x variable 
is not defined in its natural support, i.e. it can be larger than 1. Although the support 
violation is small in most models, this problem has to be considered a serious drawback 
of all CQM calculations of parton distributions, in particular if pQCD evolution of the 
model prediction is performed. We stress that our definition of the GPD H(x, £, A 2 ) can 
be easily generalized to other GPDs, which can be obtained in any CQM. 



3. RESULTS IN THE ISGUR AND KARL MODEL 



We consider the Isgur-Karl model |10| , with a proton wave function given by a harmonic 
oscillator potential including contributions up to the 2hu shell. In this case the proton 
state is given by the following admixture of states 

\N) = a s \ 2 S 1/2 )s + as'\ 2 S[ /2 )s + a M \ 2 Si/ 2 )M + <>d| 4 -Di/2)m , (7) 

where we have used the spectroscopic notation \ 2S+1 Xj) t , with t = A,M,S being the 
symmetry type. The coefficients were determined by spectroscopic properties to be: as = 
0.931, a s > = -0.274, a M = -0.233, a v = -0.067. 

The obtained behavior of the FF is shown in Fig. 1. As it is well known, such a FF 
underestimates the data at high t. Results evaluated using the IK model, Eq. (7), in the 
general formula, Eq. (4), are shown in Figs. 2 to 4. In Fig. 2, we show the t dependence 
of our results. The full line corresponds to the usual PD. One immediately realizes that 
a strong t dependence is found, in comparison with other estimates, for example, with 
the one in 0]. This has to do with the too a strong t dependence of the FF in the IK 
model. In Figs. 3 and 4 we have the full t and £ dependences. These are similar to the 
ones obtained in ||, although the £ dependence is stronger. 
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Figure 3. The GPD H(x,£,A 2 ) at A 2 = 
-0.2 GeV 2 and f = (full line), £ = 0.1 



Figure 4. The GPD H(x,£,A 2 ) at A 2 = 
-0.5 GeV 2 and f = (full line), f = 0.1 



(dotted line), £ = 0.2 (dashed line). Notice (dotted line), £ = 0.2 (dashed line). Notice 



that H(x,£, A ) is shown for x > £/2. 



that H(x,£, A ) is shown for x > £/2. 



Our results for H(x, £, A 2 ) correspond to the low momentum scale associated with the 
model. In order to compare them with the data which are going to be taken in future 
experiments, one has to evolve them to the experimental high-momentum scale. 

The proposed approach can have many interesting developments, such as the use of 
more realistic models, the inclusion of QCD evolution from the scale of the model to the 
experimental one, the addition of corrections due to relativistic dynamics, or the ones due 
to a possible finite size and complex structure of the constituent quarks, as proposed by 



several authors 11 
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